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1. Introduction 

Strange particles and otlier iieavy flavors play always a special role in the analysis of 
hadronic collisions since they carry a new quantum number not present in the incoming 
nucleons or nuclei. Therefore they are considered as one of the most promising tools 
for learning more about the dynamics of heavy ion collisions. Especially, their use as a 
signature for Quark-Gluon Plasma (QGP) formation was proposed long ago |jl|, ||. The 
argument is based on the different time scale for equilibration of strangeness due to the 
reduced kinematic threshold in the QGP [Q. There is a clear enhancement of strange 
particle production in heavy ion collisions compared to nucleon-nucleon collisions |^. 
Especially, the multistrange baryon yields increase drastically from p + p to A+A 0. 

This interpretation of the measured abundances in favor of a QGP formation 
has basically two strategies. First, dynamical arguments show that the lifetime of a 
QGP is enough to reproduce the (still not complete) strangeness equilibration [|^, |^. 
Second, the strange particle ratios are compatible with a sudden disintegrating QGP 



p!0| . However, the interpretation is still controversial and alternative explanations 



without QGP formation are also successful [p!2|, O, M. 



The various microscopic production models for strange particles in heavy ion 
collisions go from perturbative calculations on parton level |^ to string-string 
interactions (color ropes) [|T2|| to hadronic rescattering mechanisms We leave out 



the complicated question of the dynamics of strange quark production and address 
the simpler question of whether strangeness production already saturates in heavy ion 
collisions, i.e. it reaches chemical equilibrium in the flnal state. This investigation 
concentrates only on the flnal freeze-out and it implicitly assumes that the particle 
production may be described in a statistical or even thermal model. Therefore one has 
flrst to address the more general question whether the overall particle production is 
given by a thermal production mechanism. The answers will in general depend on the 
particle species, the center of mass energy and the volume. 

The recent success in the statistical interpretation of particle production in 
elementary reactions like e"*" + e~ and p + p [p + p) |]l6l as well as in heavy 



ion reactions at various energies ||17| , p!8| triggered a revival of the thermal model 
applications. Therefore we will concentrate to review the status of the thermal model 
for particle production in general and in the special case of strangeness. In section 2 
we will present a hitch hikers guide through thermal models addressing various points 
which differ in the application of the model by various groups. Section 3 is devoted 
to the discussion of our own calculation done for CERN-SPS energies. In section 4 we 
discuss the implications from the thermal analysis of particle production and give a 
general conclusion. 
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2. The thermal model 



A lot of publications have addressed the question of particle production within the 
thermal model. This can be seen in table ^ where only recent ones have been included. 
The list is certainly not complete due to ignorance or due to series of publications of one 
group where we took only the most recent one. Even if the thermal formalism is used 
for the same collision system the results vary (see table P because the model is applied 
in various approximations and extensions which we would like to discuss first. This will 
help to understand the differences in the calculations and the conclusions about the 
physics of these reactions. 



2.1. Basic particle yields 



In the thermal model the particle yields are given by a temperature T and a volume 
V common for all particles. If one considers only particle ratios then in most of the 
applications the ratio is independent of V as it was shown in p6i. In addition the 



abundance of a particle depends on its conserved quantum numbers. This is either 
regulated by chemical potentials in the grand canonical description or by restricting the 
partition function only to states which have the same quantum number as the fireball, 
i.e. canonical treatment. The basic expression for the abundance Nj of a particle of 
species j is given by 

ainZ(T,y,...) 



(1) 



All models fulfill this basic requirement with the exception of the work of Hoang for e"*" 
+ e~ collisions Some empirical formula for the particle yields is used, which is 

badly justified. Therefore all results in [11^, EO] should be taken with care. 



2.2. Statistic 

It is obvious that one should use quantum statistics for calculating the partition function, 
i.e. to use Bose and Fermi statistics. However, for practical reasons one usually switches 
to the Boltzmann approximation. The error for pions is at T = 150 MeV (T = 200 
MeV) 9.4% (11.3%), respectively and for kaons at the same temperature 0.5% (1.1%), 
respectively. This estimate suggest to use Bose statistic for pions while for all other 
heavier particles the Boltzmann approximation is valid. Note, that for entropy and 
pressure the differences between Boltzmann statistic and Bose/Fermi statistic are larger. 
Going to very low temperatures and/or high densities the use of the right statistic is 



unavoidable 37 



Table 1. Summary of calculations addressing the chemical freeze-out in various 
collisions. 
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2.3. Canonical vs grand canonical 

The strong interaction conserves exactly the quantum numbers charge Q, baryon number 
B and strangeness S. This has to be taken care in a statistical approach and therefore 
the question arises which statistical ensemble concerning these quantum numbers one 



has to use Bq, 39, 28 1. As a first estimate one usually considers the fiuctuations AO of 



a conserved quantity O in the grand canonical treatment |^ 



AO = ^(02) - (0)2 oc VN, (2) 

where is the number of all particles carrying a non zero quantum number O. In 
order to be better than 10% this suggest to use canonical treatment when the number of 
corresponding particles is below 100. However, in practice - depending on the observable 
calculated - already a much smaller total amount of particles is enough to have 
accuracy better then 10% using the grand canonical ansatz. This was shown in 



and also recently in ||2^, where particle ratios already seem to saturate for number of 



participating baryons of around 30-40. 

In the applications for heavy ion collisions the grand canonical treatment is justified 



when single particle yields are addressed |41|. We tested the difference in the resulting 



thermal parameters using canonical and grand canonical treatment for strangeness. We 
saw only minor differences between both calculations for ultra-relativistic heavy ion 
collisions. However, for strange particle correlations the use of the canonical ansatz is 
recommended IHTI . 



If the thermal model is applied to inelastic two body collisions at high energies the 
analysis should be done using the canonical ensemble. Among the calculations in table |1] 
this was only performed by Becattini |2T|, |TB[ . All other calculations concerning p + p and 
e"*" + e~ have to be taken with care. For example, we made a least mean square fit {x^- 
fit) using the grand canonical ensemble for charge, baryon number and strangeness for 
particle production in p+p collisions at y/s = 27.5 GeV. The experimental input data are 



the same as in the work of Becattini et al [|16[ and the resonances included are very much 



the same. The fit to the data gets worse but the fit temperature stays very much the 
same and turned out to be T = 160 MeV but the strangeness suppression 7^ (see Section 
|2.8|) is reduced to 7^ = 0.35 compared to 7^ = 0.51 in [|l^. The canonical treatment 
suppresses the strangeness production due to associate pair production leading naturally 
to lower strangeness yields as compared to elementary reactions. 

The chemical potentials in the grand canonical approach are first of all Lagrangian 
multipliers for the conserved charge. However, the strange quark chemical potential n., 
plays a special role in the interpretation of strange particle abundances |Q . It is zero 
in a strangeness neutral QGP but has usually non-zero values in an equilibrated hadron 
gas. The /is = line in the T-/iB plane calculated for an equilibrated hadron gas is close 
to the expected QGP phase transition line [ER EBI. This similarity leads to speculations 



s 
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about the meaning of the /Xg = hne which were extensively discussed in . 



2.4- Isospin 

Isospin breaking effects become important when target and projectile nuclei are large. In 
most of the thermal models isospin is neglected. We may estimate the effect of isospin 
breaking to be of order 1 — 2Z/A. This leads in S+Au collisions to a correction of 
order 10 % and for Pb+Pb of order 20%. Therefore one should take charge and baryon 
number separately into account especially for isospin sensitive quantities like ti^/ti^ 
or K^/K~. All above listed models in table |l] neglect isospin in the case of a heavy 
target. As an example one may discuss the tt^/tt" ratio which is different from one in 



the low pT-region of Au + Au and Pb + Pb collisions by chemical equilibration 
arguments. This ratio is sensitive to isospin violating weak decays as discussed in ||45|| . 
But in addition the isospin asymmetry in the colliding nuclei plays an important role 
and this was neglected in PH . 



The isospin SU(2) may be treated exactly like in \W\ or to a very good approximation 
as an U(l) x U(l). Then one usually conserves baryon number and charge or on valence 
quark level the net number of u-quarks and net number of d-quarks. In a grand canonical 
treatment it is equivalent to introduce chemical potential for baryon number and charge 
or for up and down quarks. 

2.5. Finite volume correction 

Calculating the partition function one switches from the summation over the states to 
the integral representation 

states „ „ 

J2 ^ J d'kgik) =A7i dk k'gik), (3) 

k 

where g{k) is the density of state. If the volume is small g{k) is very different from the 
infinite volume limit V/{2ttY and often approximated by [|^, |40[] 

V _ S L 

where V is the volume S the surface and L the circumference of the box. This is only 
an approximate formula derived for Dirichlet boundary conditions and the error is of 



the order of the last term. This kind of correction was applied in |T7, 18, 23, 29 1. 

Going to p + p collisions one might think that this kind of correction might be very 
important. However, we like to argue that in this case and also for ultra-relativistic 
heavy ion collisions it is reasonable to use the continuous density of state for the following 
reason. In equation (|^) it is assumed that the states have to fulfill Dirichlet boundary 
conditions on the edge of the interaction volume. This is only true for an infinite high 
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potential well. In reality there is no such large binding force for the particles squeezing 
its wave function to the reaction volume. The particle energies are much higher than 
a nuclear binding potential. Therefore we think that the use of the continuous density 
of state is appropriate as long as the thermal energies are above a mean field potential 
which serves as a confining box. 



2.6. Resonances 

It is now commonly agreed that the particle production is dominated by resonance 
production mechanisms. Therefore the resonance states are included in the partition 
function of all calculations in table except for |ll9| , pO[] . On the other hand their 



decay is sometimes neglected p6|, |3^ when calculations are compared to experimental 
data. Some of the publications restrict their discussion to strange baryon ratios. Then 
the resonance contribution may be restricted to the inclusion of the S*^ decay which 
feeds into the A yield. This gives a reasonable estimate for the fugacities. If, however, 
kaons are included in the chemical analysis one has to consider the sizeable feeding from 
resonances and its omission is hardly comprehensible. In the treatment of resonances 
the following items are important. 



2.6.1. Number of resonance states The modeling of a steady state hadron gas needs 
the input of all resonance states leading to the bootstrap description of thermodynamics 
of strong interacting matter . As a result one gets an exponential increasing density 



of resonance states 

p{m) (X m"' exp{m/Tii) (5) 

with the consequence of a limiting temperature Th, called Hagedorn temperature. In 
the applications for hadronic collisions it is assumed that only a finite range in the 
resonance mass spectrum is equilibrated. The finite volume of the fireball leads to a 
finite total energy giving an upper limit in the mass spectrum to be considered. This 
estimate corresponds practically to an infinite mass spectrum. It arises the question 
whether already a much lower cut in the resonance spectrum can be justified. If chemical 
equilibrium is build up via secondary interactions than the limited life time of the fireball 
leads to a limited amount of equilibrated resonance states. If, however, the hadronization 
in elementary p + p collisions follows the chemical equilibrium abundances - and this 



seemed to be the case [16| - then it is hard to argue for the omission of the higher 
mass states, especially if the temperature is around 200 MeV. Nevertheless the known 
resonance states fade at masses around 2 GeV and therefore the calculations have to be 
restricted to a finite number of states. 

The apphed cut in the resonance states of the various groups is arbitrary and given 
by practical considerations. It has been shown that there is a small influence on the 
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extracted thermal parameters on the cut in the resonance spectrum p2| , p3| , piq] . The 
temperature fitted to measured particle ratios shows a maximum at a resonance cut-off 
at ~ 1.5 GeV. It is interesting to note that fits to the experimental known density of 
states p{m) by the bootstrap formula in equation (||) start to deviate at the same mass 
of 1.5 GeV [^]. Therefore one should be aware that thermal fits at temperatures higher 
than ^170 MeV are biased from the resonance spectrum which is taken into account. 



2.6.2. Branching of resonances Not only the poor knowledge of resonance states 
around 2 GeV has some influence on the analysis for high temperatures, but also the 
branching ratios. Already at 1.5 GeV they are starting of getting basically unknown 
and the various groups use some "educated guess" . 



2.6.3. Resonance width Usually the width of a resonance is neglected. This means 
that in the Boltzmann factor exp(— v^m^ +p^/T) the mean mass m of a resonances 
is taken. This assumption was made in nearly all calculations of table |l]. For broad 
resonances this is a bad approximation. One suggestion to improve this is to distribute 



the mass states of a resonance according to a Breit-Wigner form |T5|, |T^, ^]. This 
means that the mass shell constraint in the Lorentz invariant momentum integration 



for the partition function is replaced by |49 



d> 



Pt^P 



p — 



mthr 



(6) 



where F is the width of the resonance and mthres a threshold for the resonance 
production. The inclusion of the width is important for the yield of the p-meson. It 
turned out to improve the fits in p + p collisions p6[] . 



2.7. Repulsive interaction 

The experimental particle yields and ratios are determined at (chemical) freeze-out. 
One usually assumes that the system is already such diluted that the interactions 
have effectively cease to exist. Then the problem of residual interactions don't occur. 
However, in the analysis of SPS energies the chemical freeze-out temperature seems to 
be around 160-200 MeV (see table |l|). At this temperature the particle density is still 
very high and one has to ask how density corrections influence the particle yields. In 
the bootstrap model of Hagedorn the dominant part of the attractive strong interaction 
is effectively taken into account by including the higher resonance states |^|. However, 
at high densities the repulsive part have to be accounted for, too. It is popular to use an 
excluded volume correction |ED|, pT], BB[ where the hadrons are treated as finite size 
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hard core particles. In the early suggestions |£0|, |5T| the real physical particle density 
^phy jg ];giatgcl to the ideal or point particle density n^™"* by 

nPi^y = a-VP™"^*, (7) 

where a is either given by the total point particle energy density Eq and the bag constant 
B, i.e. a = 1 + £:o/(4-B) 0], or by a; = 1 + X]j '^j'^j where is a hard core volume 



and rr^ the point particle density of species j [Q. The important point is that such 



a treatment don't influence particle ratios since the factor is common for all particle 
species and therefore don't change the thermal analysis. The volume V which appears 
as a common factor has to be regarded as the point particle volume and the physical 
volume is then given by aV . 



The above described correction is thermodynamically not consistent [50, Bl] and 



improvements have been suggested |Q Q (for an extended discussion see |^ ) . One 



may divide them basically into mean field approaches like or thermodynamically 
consistent excluded volume corrections like |5^, |53| or both . These models contain 
additional parameters which characterize either the hard core size or the mean field 
coupling Kj of a particle j. If or Kj are different for various j then they influence 
the particle ratios. The only publications in table |l] which have such an influence on 
particle ratios due to repulsions are |30|. In |^ ^ the effect is minimal since 



the mean field coupling Kb = 680 MeV fm'^ of baryons and anti-baryons is similar to 
the one of mesons Km = 600 MeV fm^. The studies in pOl show the repulsion effect 



for the heavy baryons. The extracted fiB in |^ differs from the other results of table 
|l]. Tiwari et al take for the hard core size of a particle the MIT bag model result of 
= mj/{4:B), i.e. the hard core volume scales with its mass. Therefore the massive 
baryons are additionally suppressed by their size. This has to be compensated by a 
higher fiB- 

2.8. Off- equilibrium phenomenology 

Since the life time of a fireball created in heavy ion collisions is very short and the 
dynamics is very rapid it cannot be expected that the production of particles of all 
kind follow the equilibrium statistics. In order to study the deviations from equilibrium 
quantitatively one introduces over-saturation/suppression factors 7 which measure the 
deviations from full equilibrium. For strangeness they were first introduced by Rafelski 



1 42 1 in a phenomenological way. They were defined by 



actual density 

7 = • 1°) 

equilibrium density 

It has been shown that the thermodynamically correct way of defining such a 
parameter is to define them as fugacity 

7 = exp(/i/T), (9) 
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as it is done in a grand canonical approach. This is in accordance with the model 
of relative chemical equilibrium. This means that only a subset of particles is in 
chemical equilibrium among each other and the deviation to another set of particles 
is parametrized by 7. Examples for such applications are the strangeness suppression 7^ 
the pion chemical potential or general meson and baryon suppression factors 



With the help of suppression factors one is able to study the approach to chemical 
equilibrium. Some of the calculations in table |I| allow for such a suppression and some 
don't. Since we know fromp + p collisions [^] that strangeness is suppressed by roughly 
a factor of 2 one should always allow the suppression possibility in heavy ion reactions. 
In table |l| all calculation which don't allow strangeness suppression have 1 in the column 
for 7s. 

2.9. Restricted acceptance 



Figure 1. Pion (solid line) and proton (dashed line) rapidity distribution of a static 
fireball at a temperature of T = 150 MeV normalized to one at central rapidity. 
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If thermal equilibrium is established it is either globally or locally. In most of 
the cases as for example in S+S collisions at CERN-SPS the rapidity distributions 



11 



of baryons and pions are very different in shape This suggest that the freeze- 

out parameters vary locally. One way to demonstrate this was assuming a rapidity 
dependence of the baryon chemical potential as done in [^. However, the proper 
treatment of locally varying thermal parameters is via hydrodynamics [Q. On the 
other hand, if the differences are small, the use of one global fireball in the analysis of 
4:71 data is appropriate and the result should be understood as a global average. 

The most reliable way to analyze particle yields and ratios is using 4n integrated 
data. This avoids problems due to kinematic cuts. Addressing particle ratios in 
restricted kinematic regions needs a model for the particle spectra, too. Particle spectra 
are much more sensitive to the dynamics and one needs more assumptions and knowledge 
about the space-time evolution in longitudinal and transverse direction. Therefore we 
recommend the use of An data for the study of chemical equilibration. 

Most of the analysis in table |I| is applied to particle ratios in a restricted kinematic 
region. The reason is that most of the particle ratios are only measured in a kinematic 
window. This requires that the calculation is cut to the experimental acceptance and 
the knowledge of the particle spectra is unavoidable. 

We like to demonstrate how much results in principle may depend on assumptions 
about the longitudinal dynamics at freeze-out. The two extreme scenarios are a static 
fireball and the Bjorken boost-invariant scenario ||61|]. In a static fireball the rapidity 
distribution using Boltzmann approximation is given by 
^^static Vgm^T 



dy '''' 27r2 

1 T 



2 mcosh{y) 




(10) 



while in the Bjorken case it is 

^(y) « ^Mf^m Mm/T). (11) 

It is V the volume g the spin degeneracy and m the mass. As an example we show in 
figure m the pion and proton rapidity distributions for a static fireball given by equation 
(P^). The pion and proton distribution are normalized to one and therefore its ratio at 
midrapidity is one. In the Bjorken scenario the ratio of pions to protons is given by the 
integration over rapidity in equation ([T0|) resulting in the expression of equation ([TI|). 
In our example it would be {dN^ / dy) / {dN'^ / dy) = 0.54. Since the mass of pions and 
protons are very different the effect is most pronounced in the example. We show in 
Section ^ that in practice the differences between both scenarios are minor using the 
example of S-I-Au collisions. 

The experimental rapidity distributions are broader than given by a static source 
(PUP but not infinite broad like in the case of Bjorken scaling. Since the experimental 
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width of various rapidity spectra is very similar a thermal analysis in a restricted rapidity 



range usually assumes Bjorken scaling and uses equation (|TT]) for particle yields. 



3. Results from thermal model analysis 

The strength of the thermal model is that most of the particle ratios can be explained 
by only a few parameters. However, from table |T| one cannot see how well the thermal 
model works and where the deviations start. Therefore we show the results of one 
calculation in more detail. We perform a thermal model calculation which has the 
following characteristic 



All hadronic states up to 1.7 GeV in mass are included. 

Pions follow the Bose statistic while for all other hadrons the Boltzmann 
approximation is used. 

The resonances are populated including their width according to equation (H). The 
Breit-Wigner distribution in mass is restricted to a range of two times the width 
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• When comparing to experimental results we include the feeding of resonances and 
in the case of S+Au we also include the p^-cut of the experimental ratios. 

• For S+S and S+Ag strangeness is treated in the canonical formalism while for 
S+Au we use the grand canonical ensemble. Baryon number is regulated via a 
chemical potential. Isospin symmetry is assumed. 

• No finite size correction and no repulsive interaction is included. 

We analyze experimental particle yields in two ways. On the one hand we perform 
a x^"fit to An data of S+S and S+Ag collisions and to central particle ratios in S+Au 
collisions. The second possibility is to display the experimental particle ratios in the 
T-fiB plane and look for overlap regions of the various bands. 

In table |^ we show the result of a x^"fit to 47r data from the NA35 collaboration 
(see references in table @). In S+Ag collisions the assumed isospin symmetry is slightly 
violated. In order to estimate the size of the effect we have also done a fit using separate 
chemical potentials for up and down quarks and including the total net charge. In the 
case of S+Ag collisions the result is a slightly larger = 1.572 ± 0.053 as compared to 
A„ = 1.521 ± 0.034. This is expected by the larger amount of incoming u-quarks. 



The calculations are very similar to the one of Becattini |3T| but note that in the 
case of S+S collisions a different input set of experimental data is used. Therefore we get 
a much higher temperature for S+S while for S+Ag both calculations basically agree. 
The differences between them are small deviations in the input resonances states, their 
branching and the treatment of the 77-77' mixing. 

Looking at the result in more detail one realizes that the thermal fit is not perfect, 
especially for S+S colhsions (x^/dof = 11.6/4). The largest deviations are in the anti- 
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Table 2. Result of a fit to experimental Att data of S+S and S+Ag collisions at 
CERN-SPS. The data are all measured by the NA35 collaboration. 
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^ The experimental value is extrapolated to 47r assuming the same rapidity shape as 
the A. 



baryon yields. The high absorption cross section of this particles may explain the 
deviations. 

A different way of displaying the quality of the thermal model approach is shown 
in figure ^ where in the T-/iB-plane various bands indicate experimental particle ratios. 
We changed to particle ratios because they are nearly independent of the volume. Since 
we use the canonical ensemble for strangeness we have a small influence of the particle 
ratios on the volume. We calculated the used experimental particle ratios from table 
The errors of the ratios are determined by adding the individual errors quadratically. 
The bands correspond to the upper and lower bound on the experimental ratio. Note 
that a fixed volume was used. 

In figure |^ we see no real overlap region of all particle ratios. Especially, the ratios 
containing the h~ fail to cover the y^-^i point which is given by the filled circle. We 
like to point out the possible sign of an enhanced entropy production seen in the h~ 
as discussed in [0, Our present reevaluation of the experimental data confirms 
this possibility. We expect a stronger effect on an enhanced pion production in Pb+Pb 
collisions. 

In S+Ag collisions the quality of particle production using a thermal model is similar 
to the one in S+S as may be seen from the table |^ or in figure |^ where we plotted 
again various ratios in the same way as in figure 0. All particle ratios, excluding K^/h~ , 
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Figure 2. Experimental particle ratios in the T-^b plane for S+S collisions taking 
7s = 0.84 and V = 81.5 fm'^. The bands along the abscissa correspond to the following 

ratios (going from left to right): jh^ ( ), A/(p - p) ( ), A/A ( ), 

VI{V-V) ( ), h-/{p-p) ( ) K+/K- ( ) and K^J K ( ). The 

point indicates the result of the x^-fit in table 0. 




have one common overlap at T = 170 ± 10 MeV and = 220 ± 10 MeV. The inclusion 
of the I h~ ratio in the x^-fit moves that result out of the otherwise common overlap. 
Excluding K'^ one gets for S+Ag collisions similar freeze-out parameters as for S+Au 
collisions below. 

There are not enough An data on S+Au collisions and therefore we switch to particle 
ratios. The analysis is inspired by the work of Braun-Munzinger et al |18|. However we 



take in our analysis only a subset of particle ratios from their list, excluding all ratios 
which don't cover midrapidity ?/cm = 2.65. For the particle yields we use the scaling 
assumption, i.e. equation (|lTD as it was done in |T^. In addition we change to the grand 
canonical ensemble for strangeness. 

The result of the x^-^t is given in table ^ We reproduce the temperature T = 160 
MeV as it was assumed in |T^ but we got a slightly lower /x^ of 158 MeV. The main 
difference of our calculation to the one in [jl8| is that we allow strangeness suppression. 
The ratios sensitive to 7^ are S~/A and i?n = (^2 + + H^). The resulting low 

value of 7s = 0.65 is not in agreement with the assumption of full chemical equilibrium 



15 



Figure 3. Experimental particle ratios in the T-fiB plane for S+Ag collisions taking 
7s = 0.82 and V — 275 fm'^. The bands along the abscissa correspond to the following 

ratios (going from left to right): K'^/h" ( ), A/{p~p) ( ), A/A ( ), p/{p-p) 

( ) and /{p — p) ( — • — ). The point indicates the result of the x^"fit in table 




for strangeness as it was assumed in many calculations. 

The result of table |^ shows that a thermal hadron gas model even with no complete 
strangeness equilibration is not able to reproduce all experimental data. Some serious 
deviations are not in the table like S/A = 0.23 ± 0.02 which in the thermal model 
is given by H/A = 0.135. The disagreement in the multi-strange baryons might indicate 
the onset of non-equilibrium physics with the origin in a QGP formation |l[]. This 
proposal was recently discussed in depth in [|]. 

We tested our result against the assumption of Bjorken scaling in rapidity and did 
the same fit assuming one static fireball, i.e. using equation (^) at y = 0. We got a 
rather similar result of T = 158 ± 3 MeV, Xq = 1.408 ± 0.012 and 7, = 0.74 ± 0.05. We 
see in the x^-fit no large dependence on the assumption about the rapidity distribution. 

We show in figure ^ some selected ratios from table ^ together with the point from 
the x^-fit. A gain we see no perfect agreement but a rather broad region where various 
bands concentrate. 
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Table 3. Result of a fit to experimental central rapidity particle ratios in S+Au(W,Pb) 
collisions at CERN-SPS assuming Bjorken scaling. The extracted thermal parameters 
are T = i60.2 ± 3.5 MeV, A, = i.390 ± 0.012 and 7^ = 0.656 ± 0.04i. It is 
xVdof = 37.7/10. 



S+Au 



q 

p/tt- 
P/P 

K+/K- 

K'JA 

A/A 

n/n 



cal. 



data 



target rapidity prcut ref. 



0.0782 

0.188 

0.0262 

0.139 

0.0816 

2.03 

1.57 

1.21 

0.203 

0.0967 

0.283 

0.145 

0.430 



0.088 ±0.007 
0.19 ±0.03 
0.024 ±0.009 
0.12 ±0.02 
0.15 ±0.02 
2. 14 ±0.06 
1.67±0.15 
1.4 ±0.1 
0.196 ±0.011 
0.097 ±0.006 
0.47 ±0.06 
0.8 ±0.4 
0.57 ±0.41 



Pb 

Pb 

Pb 

Pb 

Au 

W 

W 

W 

W 

W 

W 

W 

W 



2.3-3.0 

2.6-2.8 

2.6-2.8 

2.65-2.95 

2.1-2.9 

2.5-3.0 

2.3-3.0 

2.5-3.0 

2.3-3.0 

2.3-3.0 

2.3-3.0 

2.5-3.0 

2.5-3.0 









1.0 
0.9 
1.0 
1.2 
1.2 
1.2 
1.6 




" Dg = {h+ -h-)/{h+ +h-) 

Rk = {K+ +K-)/K°^ 



EMU05 
NA44[|8) 
NA44[|| 
NA44g 
WA80 
WA85 
WA85 
WA85 
WA85 
WA85 
WA85 
WA85 
WA85 



The statistical error in the x^-fit is determined by the region where increases by 
one unit from its minimum. However, from the figures 0, ^ and ^ one sees that the 
systematic error of the model applied is much larger. From the figures we conclude that 
the freeze-out temperature in heavy-ion collisions is still very uncertain and one has in 
fact to take a range of T^^^"^ = 150-200 MeV. 

The results on sulphur induced colhsions at CERN-SPS indicate no full strangeness 
equilibration but one is very close to it. The slightly higher 7^ in the smaller collision 
systems may be a result of having no multi-strange particles in the x^"fit- 



4. Conclusions from thermal models 



Since the application of a statistical model to multi-particle production by Fermi ||7^ its 
validity is under debate. Therefore we would first like to make some general remarks. 
Thermodynamics is first of all a formalism which may be derived from statistical 
(quantum) mechanics in the infinite volume limit. The basic assumption going in is 
that all states which are allowed by conservation laws, including energy conservation, 
are equal probable. We get the microcanonical formalism. Going to the canonical 
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Figure 4. Experimental particle ratios in the T-^b plane for S+Au(W,Pb) collisions 
taking 7^ = 0.66. In this plot the experimental kinematic cuts are not accounted for. 
The bands are explained along the upper part of the figure from left to right: p/vr^" 

( ), i^o/A (_" --),D,{ ), K+/K- ( ), A/A ( ), ( ) 

and the broad VL/VL ( ). The point indicates the result of the x^"fit in table |3[ 




formalism only states sharing the same total energy have the same probability which is 
proportional to oc exp(— i?/T). The equal probability is usually violated in elementary 
reactions where the final state probability is given by the corresponding matrix element. 
Therefore the general opinion is that in elementary reaction the thermal model has no 
justification. 

However, a large number of particles in the final state suggests to use a statistical 
approach. Therefore the thermal model was applied for elementary reactions and 
one gets reasonable agreement [|1^, This observation suggest that if only enough 
energy is available the particle production is dominated by the statistical component 
and dynamical aspects are of minor importance. In addition Becattini |15, 16 1 got 
the important result that there is an universal hadronization temperature in different 
elementary collision systems and it is independent on ^/s. The interpretation is that in 
the rest frame of the leading particle/parton the probability of producing a particle with 
energy E is proportional to oc exp(— E'/T/^). The hadronization temperature Th may be 
identified with the Hagedorn temperature Th- In the interpretation of Hagedorn it 
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is not possible to create a system of higher temperature than Th unless there is a phase 
transition. The limiting temperature is seen in elementary reactions even at very high 



energies like in p+p at CERN |T6[. The above described observation was also dubbed 
as statistical filling of phase space. We want to point out that there is a difference in 
the basic hadron production probability compared to the string models for hadronic 
reactions ||31|. There the production probability goes basically like oc exp(— * k) 
with k being a universal constant. 

The mechanism for chemical equilibration in nucleus-nucleus collisions is expected 
to be different from the one in elementary reactions. In nuclear reactions we assume 
that chemical equilibration is established by secondary interactions among the produced 
hadrons. Therefore we distinguish two mechanisms which bring the system to maximum 
entropy or chemical equilibrium. 

• The production of particles, i.e. the hadronization or fragmentation, follows a 
statistical law and already at their production they are distributed according to 
maximum entropy. The ensemble average is done by averaging over many events 
in the experimental analysis. 

• The maximum of entropy is build up in the classical sense by interactions among 
the particles until detailed balance is reached. The ensemble average is reached in 
each collision by the average over the lifetime of the system (engodic theorem). 

The thermodynamical formalism cannot distinguish between both scenarios and one 
has to use dynamical arguments to justify one or the other mechanism. 

Since at SPS-energies the chemical freeze-out temperatures and the chemical 
potentials are very similar in p+p collisions and in S+A collisions (see table |1|) one cannot 
use them to justify secondary interactions for chemical equilibrium. A superposition 
of p+p collisions explains most of the features in the nuclear collisions [76| with the 
exception of strangeness. Therefore we emphasize the importance of the measurement 
of strangeness because there the difference in p+p to A+A is most clearly seen. 

The freeze-out temperature is expected to decrease with increasing A, because 
freeze-out occurs when the mean free path is of the order of the size of the system. Such 
an effect is not observed in an unambiguous way, yet. So far only indications are seen 
as for example the decrease of chemical freeze-out in our analysis from S+S to S+Ag to 
S+Au. A clear sign for chemical equilibration due to secondary interactions would be a 
difference in the chemical freeze-out of p+p compared to a real heavy nucleus like Pb+Pb 
or Au+Au. Such an analysis hasn't been performed yet but with the now analyzed data 
of Pb+Pb at CERN-SPS (see for example the various contributions to this proceedings) 
it will be possible soon. At the AGS we expect different freeze-out temperatures for 
Si+Au and Au+Au. The first results on Au+Au go in this direction but we have 
to wait for the completion of the experimental data analysis. 

Coming finally back to strangeness we have here a clear signal of the difference 
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between elementary reactions and nuclear collisions. The strangeness suppression factor 
7s is the quantitative measure for the strangeness enhancement in the thermal model. 
In p+p it is 7s ~ 0.5 [|T^ but in nuclear collisions at CERN-SPS it is around 0.7-1 (see 
table |I|). A strangeness enhancement is seen at the AGS, too |7^, but the situation is 
not so clear in terms of the thermal ansatz. A consistent study of the 7^ dependence 
has not been made. First there is no thermal analysis of the p+p interaction at the 
corresponding y/s and second the analysis at AGS assumes mostly full strangeness 
equilibration as it is seen in table |I]. However, we have already remarked [rT| that a 
7s ~ 0.7 is better for describing the data at AGS. 

We have shown that the thermal model fits are not perfect including all measured 
particle species. The deviations are a source of debate. The various interpretations are 
that the thermal description is not valid at all, the hadronization of a QGP leaves non- 
equilibrium tracks in special hadron ratios [H, anti-baryons exhibit a large absorption 



|78| or the deviations are not serious |]18[. The final answer will be given in the expected 
high statistic data of Pb+Pb in the future. 

We summarize that there are strong signs of chemical equilibration in heavy ion 
collisions. Since there is already an equally good chemical equilibrium (excluding 
strangeness) in the basic p+p collisions chemical equilibrium cannot be used for 
justifying secondary hadronic collisions. (There are better signals for abundant 
secondary interactions like the collective flow studies fT^.) However, the strangeness 
production is very different between both collision systems and it may be used as the 
chemometer for chemical equilibrium. 
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Additional Figure 



Figure 5. Overview of table |l| in the r-/iB-plane. The plot uses color coding and 
information is lost using black and white. The plot is not included in the version 
appearing in Journal of Physics G 
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